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A small body moving in the field of a much larger black hole and subjected to its own gravity moves 
on an accelerated world line in the background spacetime of the large black hole. The acceleration 
is produced by the body's gravitational self-force, which is constructed from the body's retarded 
gravitational field. The adiabatic approximation to the gravitational self-force is obtained instead 
from the half-retarded minus half-advanced field. It is much easier to compute, and it is known to 
produce the same dissipative effects as the true self-force. We argue that the adiabatic approximation 
is limited, because it discards important conservative terms which lead to the secular evolution of 
some orbital elements. We argue further that this secular evolution has measurable consequences; 
in particular, it affects the phasing of the orbit and the phasing of the associated gravitational wave. 
Our argument rests on a simple toy model involving a point electric charge moving slowly in the weak 
gravitational field of a central mass; the charge is also subjected to its electromagnetic self-force. 
In this simple context the true self-force is known explicitly and it can cleanly be separated into 
conservative and radiation-reaction pieces. Its long-term effect on the particle's orbital elements can 
be fully determined. In this model we observe that the conservative part of the self-force produces 
a secular regression of the orbit's periapsis. We explain how the conclusions reached on the basis of 
the toy model can be extended to the gravitational self-force, and we attempt to extend them also 
to the case of rapid motions and strong fields. While the limitations of the adiabatic approximation 
are quite severe in a post-Newtonian context in which the motion is slow and the gravitational field 
weak, they may be less severe for rapid motions and strong fields. 



I. INTRODUCTION 

The gravitational inspiral of a solar-mass compact ob- 
ject into a massive black hole residing in a galactic center 
has been identified as one of the most promising sources 
of gravitational waves for the Laser Interferometer Space 
Antenna . The need for accurate theoretical models of 
the expected signal, for the purposes of signal detection 
and source identification, has motivated an intense effort 
from many workers to determine the motion of the small 
body in the field of the large black hole. This is done in a 
treatment that goes beyond the geodesic approximation 
and takes into account the body's own gravitational field, 
which is a small perturbation over the field of the black 
hole; this must be done without relying on slow-motion or 
weak-field approximations. In this treatment the small 
body can be thought of as moving on a geodesic of a 
perturbed spacetime, or equivalently, it can be thought 
of as moving on an accelerated world line in the back- 
ground spacetime of the large black hole. The second 
point of view is generally adopted, and the body is said 
to move under the influence of its own gravitational self- 
force 0, Q • The self-force is derived from the retarded 
gravitational perturbation produced by the moving body. 
For a review of the self- force formalism, see Ref. Q|. 

The concrete evaluation of the gravitational self-force 
acting on a small body moving in the Kerr spacetime 
is a challenging project that has not yet been com- 
pleted. Part of this challenge is concerned with the re- 
construction of the metric perturbation 0, 0> from 
the Teukolsky variables, which can be more practically 
evaluated [8(. Another is concerned with the regular- 
ization of the body's retarded field near the world line 



0, 0, , in a context where it is expressed as 

an infinite sum over spherical- harmonic modes. Yet an- 
other challenge resides in the fact that the gravitational 
self- force is a gauge-dependent quantity ; this implies 
that the improved equations of motion must be incorpo- 
rated in a self-consistent wave-generation formalism be- 
fore they can yield gauge-invariant waveforms. 

II. ADIABATIC APPROXIMATION. . . 

Given this state of affairs, it is tempting to seek ap- 
proximations to the self-force formalism that may by- 
pass some of these challenges but still produce accept- 
ably accurate results. One such approximation was for- 
mulated by Mino 0EII3, who showed that the long- 
term evolution of the three principal orbital elements (the 
body's orbital energy, angular momentum, and Carter 
constant, which would all be constant under geodesic mo- 
tion but instead evolve under self-forced motion) can be 
reproduced on the basis of a radiation-reaction force con- 
structed from the half-retarded minus half-advanced so- 
lution to the perturbation equations. In particular, Mino 
showed that the total work done by the radiation-reaction 
force equals the total gravitational-wave energy radiated 
by the moving body, a result that was previously estab- 
lished by Quinn and Wald [II] and Gal'tsov Mino 
was further able to show that these long-term evolutions 
are, to a very good approximation, gauge- invariant. 

Because the retarded and advanced fields are equally 
singular on the world line, their subtraction produces 
a regular field, and the computation of the radiation- 
reaction force requires no regularization. This compu- 
tation is therefore much simpler than the calculation of 
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the true self-force, and the gauge-invariant nature of the 
radiation-reaction force argues in favor of its direct in- 
volvement in waveform calculations. 

The approximation of the true self-force by the 
radiation-reaction force discards the rapid, oscillatory 
changes of the three principal orbital elements which oc- 
cur on the time scale of the orbital period; it retains 
only the slow changes that survive after time-averaging. 
This replacement of the true self-force by the radiation- 
reaction force in a calculation of the orbital evolution is 
called the adiabatic approximation. The computation of 
inspirals in the adiabatic approximation, and of their as- 
sociated gravitational-wave signals, has so far been pur- 
sued by at least two research groups [2(], 0, 03 . 

Confidence in the reliability of the adiabatic approxi- 
mation rests on the results of Mino (which are fully re- 
liable) and the belief that all secular effects associated 
with the true self-force are embodied in the long-term 
evolution of the three principal orbital elements (energy, 
angular momentum, and Carter constant). Because the 
true self-force differs from the radiation-reaction force 
by purely conservative terms (which do no work), this 
is a statement of belief that the conservative part of 
the true self-force produces only short-term effects that 
do not accumulate over time and do not survive after 
time-averaging. This belief is expressed, for example, by 
Drasco, Flanagan, and Hughes |20j, who state near the 
beginning of their Sec. 1.2: "The effect of the dissipative 
pieces of the self force will accumulate secularly, while 
the effect of the conservative pieces will not. Hence the 
effect of the dissipative pieces on the phase of the orbit 
will be larger than that of the conservative pieces by a 
factor of the number of cycles of inspiral." 



III. . . . AND ITS LIMITATIONS 

Our purpose in this paper is to show that this belief is 
not well founded: The conservative terms in the true self- 
force do produce effects that accumulate over time and 
significantly affect the phasing of the orbit. Our meth- 
ods allow us to conclude that these effects are important 
in slow-motion, weak-field situations. They do not, how- 
ever, allow us to investigate directly more interesting sit- 
uations involving rapid motions and strong fields. Nev- 
ertheless, a tentative extrapolation of our slow-motion 
results indicate that these effects may be less important 
in strong-field situations. 

While the radiation-reaction force captures the secu- 
lar evolution of the three principal orbital elements (en- 
ergy, angular momentum, and Carter constant), it makes 
no statement regarding the evolution of the remaining 
three orbital elements. (We shall call these, for reasons 
that will become clear presently, the positional orbital el- 
ements.) In particular, the radiation-reaction force does 
not capture the secular evolution of the positional ele- 
ments. 

In the case of geodesic motion all six orbital elements 



are constants of the motion. The three principal elements 
serve to label the geodesic, and the three positional ele- 
ments serve to specify the body's initial position on this 
geodesic. Under orbital evolution driven by the true self- 
force, all six elements vary with time, and the motion is 
tangent to an evolving osculating geodesic. While it is 
clear that the principal elements must evolve secularly, 
we intend to show that the positional elements also can 
evolve secularly, and that this evolution is driven by the 
conservative part of the true self-force. The adiabatic 
approximation, therefore, does not account for the secu- 
lar evolution of the positional elements, and we intend to 
show that this leads to measurable phasing effects in the 
waves. 



IV. ELECTROMAGNETIC SELF-FORCE IN A 
WEAKLY CURVED SPACETIME 

Our argument is based on a simple toy problem in 
which the true self-force can be evaluated explicitly 
and decomposed cleanly into conservative and radiation- 
reaction pieces. The toy problem involves a point electric 
charge q of mass m moving slowly in the weak gravita- 
tional field of a (noncompact) star of mass M. In the 
toy problem the charge is a substitute for the small com- 
pact body, the electromagnetic radiation is a substitute 
for the gravitational radiation, and the star is a substi- 
tute for the massive black hole. We take the star to 
be immobile and we place it at the origin of a three- 
dimensional Cartesian coordinate system x = [x, y, z]. 
The electromagnetic self-force acting on the point charge 
was calculated by DeWitt and DeWitt [24[, building on 
the foundations laid by DeWitt and Brehme (2{| ; this cal- 
culation was reviewed recently in Ref. [2(|. Adopting a 
three-dimensional, Newtonian language throughout this 
section, the equations of motion are 

d 2 x 

m -Jf2 = m 9 + /self, (1) 

where x(t) is the charge's position vector, 



is the Newtonian gravity of the central mass (with r = \x\ 
and r = x/r) , and 

/self = Ac^r + Kr\q 2< ^ (A c = A rr = 1) (3) 

is the electromagnetic self-force. The first term on the 
right-hand side of Eq. is the conservative part of the 
self-force; the second term is the usual expression for the 
electromagnetic radiation-reaction force. We work with 
units such that G = c = 1, and we have inserted parame- 
ters A c = 1, A rr = 1 in order to later distinguish between 
conservative and radiation-reaction effects. In Appendix 
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A we show that Mino's prescription [15( for the radiation- 
reaction force gives rise precisely to the second term of 
Eq. • The adiabatic approximation to the electromag- 
netic self-force therefore consists of setting A c = and 
A rr = 1 in Eq. ©. 

Under the action of g alone the point charge would 
trace a Keplerian orbit of semilatus rectum p and eccen- 
tricity e, and its orbital radius would be described by 



P 



1 + e cos(0 — u>) ' 



(4) 



where (f> is the charge's true longitude and the constant 
lu is the longitude at periapsis. In Eq. @ only r and 
4> depend on time. The charge's energy per unit mass 
is E = — M (1 — e 2 )/(2p) and its angular momentum per 
unit mass is L = ^/Mp. The orbital elements (p, e) there- 
fore act as a substitute for (E, L) and can be adopted as 
the principal orbital elements. The remaining element u> 
is a positional orbital element. [A second positional ele- 
ment, to, does not appear explicitly in Eq. Q and will 
not be needed in the sequel; it is defined by the state- 
ment 4>(to) = uj.] Because the self- force of Eq. (01 keeps 
the orbital motion in a fixed plane, there is no need to 
introduce additional orbital elements. 

Under the perturbation produced by the self-force the 
orbital motion is still described by Eq. |@J but the ele- 
ments (p, e, uj) acquire a time dependence; the motion is 
tangent to an evolving osculating Keplerian orbit. Em- 
ploying standard methods from celestial mechanics we 
show in Appendix B that the long-term evolution of the 
orbital parameters is governed by the differential equa- 
tions 



(p) = -A rr 



4q 2 M (1 - e 2 ) 



2\3/2 



2\3/2 



= -A, 



P 

q 2 M e(l - e 

m p a 

lq 2 M f~p~ ( 1 - e 2 f' 2 
p3 



(5) 
(6) 
(7) 



Here the overdot indicates a time derivative, and the an- 
gular brackets mean that (p, e, to) are averaged over a 
complete period of the unperturbed orbit. It is assumed 
that q 2 <C mp throughout the evolution, so that the time 
scales associated with the changes in (p, e, uj) are much 
longer than the orbital period. 

Equations 10- describe the secular evolution of the 
orbital elements; they average out the unimportant os- 
cillatory behavior that leads to a zero cumulative change 
over a large number of orbits. The decay of the princi- 
pal elements p and e describes a shrinkage of the orbit 
accompanied by circularization. The presence of A rr = 1 
in Eqs. JSJ and I© shows that these effects are purely 
dissipative, and that the evolution of p and e is driven 
entirely by the radiation-reaction part of the self-force. 
This evolution would be reproduced exactly by the adi- 
abatic approximation. The negative sign on the right- 
hand side of Eq. indicates that the orbit undergoes 




time : initial period) 

FIG. 1: Magnetic field produced by the orbiting particle, 
evaluated as a function of time in the wave zone, within the 
electric-dipole approximation. The units of the magnetic field 
are arbitrary, and the time is given in units of the initial 
orbital period. The initial orbital elements are po — 20M, 
eo = 0.5, and ujq — 0. We set q 2 /(mpo) = 0.002. In the color 
version of the figure, the blue curves represent the wave gener- 
ated by an orbital evolution driven by the true self-force, and 
the red curves represent the wave calculated within the adia- 
batic approximation. In the black-and-white version the true 
wave appears thick and dark, while the approximate wave ap- 
pears thin and light. The main plot shows the entire inspiral. 
The fact that the red/light curve appears to move upward 
relative to the blue/dark curve is a consequence of a relative 
phase shift driven by the secular evolution of uj. The lower 
inset on the left shows the earliest part of the inspiral, before 
the phase shift had a chance to accumulate. The lower inset 
on the right shows a later portion of the inspiral, with the 
approximated wave (red/light) leading slightly in phase with 
respect to the true wave (blue/dark). Finally, the upper inset 
shows the latest portion of the inspiral, when the two waves 
are very much out of phase. 



also a regression of its periapsis; the presence of A c = 1 
shows that this effect is purely conservative. The evolu- 
tion of uj is driven entirely by the conservative part of the 
self-force, and this evolution would not be reproduced by 
the adiabatic approximation. 

We see that the conservative term in the electromag- 
netic self- force is indeed responsible for a secular effect, 
the accumulating regression of the orbit's periapsis. This 
effect would be missed altogether by the adiabatic ap- 
proximation, which fails to account for the secular evo- 
lution of the positional orbital elements. 

The secular evolution of uj has a direct effect upon the 
phasing of the orbit, and upon the phasing of the as- 
sociated electromagnetic-wave signal. This statement is 
illustrated in Fig. 1, which displays the electromagnetic 
wave generated by the orbiting particle; the calculations 
behind this figure are presented in Appendix C. The fig- 
ure shows very clearly that a model waveform calculated 
on the basis of the adiabatic approximation (which ig- 
nores the evolution of uj) gradually goes out of phase 
with the true signal. This is in spite of the fact that both 
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waves incorporate the same correct long-term evolution 
of p and e. 

The time scale associated with the dephasing can be 
estimated on the basis of Eq. J2J) . Assuming that e is not 
too close to unity, this is given by (uj) ~ 1/Td p h, or 



Tdph 



M 



5/2 



(8) 



The radiation-reaction time scale, on the other hand, can 
be estimated on the basis of Eq. 10. This is given by 
(p) ~P/r IT , or 



(9) 



The dephasing time is shorter than the radiation- reaction 
time by a factor of order y/ M/p <C 1. The total phase 
shift A$ accumulated during a radiation-reaction time is 
estimated as A<£> ~ (u})t ti ~ r rr /rd p h, or 



A$ ~ ^/p/M > 1 



(10) 



after using Eqs. (jHJ) and JHJ). The accumulated phase shift 
is large, and this indicates that the adiabatic approach 
will not produce an adequate approximation of the true 
waveform over such a time interval. 



V. DISCUSSION 

Our two main conclusions are these: 

• The secular evolution of the orbital elements is not 
driven only by the radiation-reaction part of the 
self-force; the conservative terms participate also. 
In particular, the conservative part of the self-force 
gives rise to a secular evolution of the positional 
elements. 

• The secular evolution of the positional elements, 
which is not accounted for by an adiabatic ap- 
proximation to the self-force, has measurable con- 
sequences. In particular, this evolution affects the 
phasing of the orbit and therefore the phasing of 
the associated wave. 

These conclusions apply to the toy problem presented in 
the preceding section, but we firmly believe that they are 
not limited to this specific example. We are convinced 
that they hold in the context of the gravitational self- 
force acting on a small body moving in the field of a 
Kerr black hole. 

Returning to the gravitational problem, we emphasize 
that it should not come as a surprise that the conservative 
terms in the self-force can drive the secular evolution of 
the positional orbital elements. And it should not come 
as a surprise that this evolution can produce a measure- 
able effect on the phasing of the gravitational waves. We 



shall elaborate on this observation in the next two para- 
graphs, but we note first that a similar point was made 
recently by Ajith et al. [H|l|. 

Consider the post-Newtonian description of a two- 
body system, in a regime where m (the small mass) is 
much smaller than M (the large mass) . The equations of 
motion for the relative orbit take the schematic form 



a — <70PN + <7lPN + <72PN + 92.5PN 



(11) 



in which a is the acceleration of the relative position vec- 
tor, <7opn stands for the Newtonian gravitational field, 
and the additional terms are labeled by their post- 
Newtonian order. Each term in the post-Newtonian ex- 
pansion can be expanded in powers of m/M, so that 



o 1 

9nPN = SnPN + Jj9 n PN 



0(m 2 ). 



(12) 



Summation over post-Newtonian orders of the terms 
<7nPN gives rise to an acceleration a which describes the 
geodesic motion of a test mass in the gravitational field of 
the mass M. Summation of the terms (m/M)g^ PN gives 
rise to (m/M)a 1 , the correction to the geodesic motion 
that comes from the true self-force. These considerations 
show that in a post-Newtonian context, the self- force has 
conservative terms at order OPN, 1PN, 2PN, and higher, 
and it has dissipative terms at order 2.5PN and higher. 

It is now easy to see why the conservative part of the 
self-force must drive a secular evolution of the positional 
elements. It is known, for example, that the test-mass 
term gJpN drives a secular shift in u> relative to Keple- 
rian motion; this is exemplified by the famous perihelion 
advance of Mercury. It is also known that (m/M)gj PN , 
the self-force term, contributes a correction of order m/M 
to this shift |23 ■ Concretely, the post-Newtonian expres- 
sion for the periapsis advance is given by Eq. (12.29) of 

Ref. laH, 



(W/1PN 



3(1 -e 2 ) 3 / 2 (M + m) 3 / 2 



P 



5/2 



(13) 



since this depends on the total mass M + m, there is a 
test-mass contribution at order ra° and a self-force con- 
tribution at order m. The self-force therefore produces 
a secular shift in u), over and above the effect already 
present in the test-mass limit. Because the conservative 
part of the self-force begins at a low post-Newtonian or- 
der, it should be clear that it will drive an evolution 
that is potentially more significant than the evolution 
driven by dissipative effects, which begin at order 2.5PN. 
This property was featured in the toy model: The con- 
servative part of the electromagnetic self-force can for- 
mally be thought of as a 1PN term in the equations of 
motion, while the radiation-reaction force is a smaller 
1.5PN term. As we saw this led to a dephasing time that 
was shorter (by a factor of the orbital velocity) than the 
radiation-reaction time. 

Incidentally, an extension of this argument shows why 
a recent claim made by Burko |3lj | must be erroneous. 
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Burko claims that if the mass m is spinning rapidly, then 
the effects of the conservative spin-orbit interaction will 
dominate over conservative self-force effects. It is easy 
to see why this claim must be false: Formally the spin- 
orbit contribution to Eq. Ijlll) begins at 1PN order, but 
for rapidly-spinning compact bodies the pos t-Newtonian 
scaling is in fact shifted to 1.5PN [32L l33|. Thus the 
spin-orbit interaction contributes a term gi.spjsi, and the 
fact that the spin angular momentum scales as the mass 
squared implies that g sp i n -orbit = (m/M )g\ 5PN . This 
addition to the self-force is therefore of 1.5PN order, and 
it will not dominate over the 1PN term that produces 
the periapsis shift. 

These observations generalize to rapid motions and 
strong gravitational fields. In this context the conserva- 
tive part of the self- force should still be expected to drive 
a secular evolution of the positional orbital elements — 
those that characterize the initial position of the small 
body on the osculating geodesic. The secular evolution 
of the positional elements will affect the phasing of the or- 
bit and the phasing of the associated gravitational wave. 
It is not, however, accounted for by the adiabatic ap- 
proximation, which captures only the dissipative aspects 
of the orbital evolution. 

We must add that there are important differences be- 
tween weak-field and strong-field situations, and admit 
that our weak-field methods do not allow us to determine 
the strong-field effects associated with the conservative 
part of the self-force. We shall, nevertheless, attempt 
to draw some conclusions by extrapolating our weak- 
field results to strong-field situations. These conclusions, 
of course, are tentative and await confirmation from a 
proper strong-field investigation. 

Our first observation is that in the case of rapid mo- 
tions and strong fields, the conservative-driven evolution 
occurs over a time scale that is now comparable with 
(and not much shorter than) the radiation-reaction time 
scale. The separation of scales no longer occurs because 
the post-Newtonian scaling with the orbital velocity does 
not apply when the velocity is close to unity. This im- 
plies that the total accumulated phase shift now amounts 
to a number of wave cycles that may not be much larger 
than unity. Indeed, the estimate of Eq. (|l(Jf> must be 
replaced by A<I> ~ (p/M ) 3 / 2 in the case of gravity; this 
number would be large in a post-Newtonian situation, 
but in the present context we have p/M > 1 and this 
leads to A<E> > 1. This represents a fractional correc- 
tion of order m/M to the total number of wave cycles 
accumulated during the entire inspiral, which is of order 
M/m. 

This strong-field estimate of the size of the accumu- 
lated phase shift is in agreement with statements made in 
Ref. [2l| and in Sec. 1.2 of Drasco, Flanagan, and Hughes 
po| . Recall the statement quoted earlier in Sec. II: "The 
effect of the dissipative pieces of the self force will ac- 
cumulate secularly, while the effect of the conservative 
pieces will not. Hence the effect of the dissipative pieces 
on the phase of the orbit will be larger than that of the 



conservative pieces by a factor of the number of cycles of 
inspiral." While we still disagree with the first sentence, 
we acknowledge that the second sentence may well be 
true in the context of strong fields and rapid motions. 
But the fact that A<E> scales as (p/M) 3 / 2 and must there- 
fore be large in post-Newtonian situations seems to have 
been overlooked by these authors. 

As a final remark we add that our considerations are 
guided by the realistic expectation that the orbital mo- 
tion of a small compact body around a massive black 
hole will be eccentric and taking place outside the equa- 
torial plane of the rotating black hole. For the special 
case of equatorial, circular orbits the relevance of the 
positional orbital elements disappears; the secular drift 
of to, for example, has no measurable consequences for 
circular orbits. In such special cases the adiabatic ap- 
proximation can be expected to be reliable. We suspect 
that it is largely the consideration of circular, equatorial 
orbits that has boosted the confidence in this approach; 
refer, for example, to Appendix A of Drasco, Flanagan, 
and Hughes pq . 
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APPENDIX A: RETARDED, ADVANCED, AND 
RADIATION-REACTION FORCES 

The electromagnetic self-force acting on a point charge 
q moving in an arbitrary curved spacetime was first cal- 
culated by DeWitt and Brehme |25l , and their expression 
was later corrected by Hobbs [34|]. We briefly sketch the 
main steps of the derivation here, relying heavily on the 
presentation given in Sec. 5.2 of Ref. [j] (hereafter re- 
ferred to as LRR). We then provide an expression for 
the advanced version of the self-force. Subtracting the 
two gives Mino's radiation-reaction force [13 . Finally, 
we evaluate these forces in the case of a charge moving 
slowly in the weak gravity of a central mass M, rely- 
ing heavily on the work of DeWitt and DeWitt [24| as 
reviewed in Ref. [2(| (hereafter referred to as PP). 

The derivation of the true, retarded, self-force begins 
with the following expression for the retarded potential 
produced by a charged particle moving on a world line de- 
scribed by the parametric relations z m (t) [LLR Sec. 5.2.2, 
Eq. (443)]: 

A« t (x) = q J G Ic ^(x,z(r))u^r)dT. (Al) 
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Here x is the spacetime point at which the potential is 
evaluated, G Te " a ,(x,x') is the retarded Green's function 
associated with the wave equation satisfied by A a (in the 
Lorenz gauge), and u M = dz^/dr is the charge's velocity 
vector. 

The electromagnetic field -F r "f (x) obtained from the re- 
tarded potential is singular on the world line, and it must 
be regularized before taking the limit x — > z and evaluat- 
ing the Lorentz force qF^ v {z)u v . Detweiler and Whiting 
|!j have shown that the correct regularization procedure 
is to remove from F r "f (x) a singular field F^(x) which 
is known to exert no force on the particle. The singular 
field is defined and evaluated in Sec. 5.2.5 of LRR, and 
the regular remainder 



Ff(x)=F I a / t (x)-Ff(x) 



(A2) 



is what must be substituted into the Lorentz-force equa- 
tion. 

When the charged particle moves freely (being sub- 
jected to no other force but its own self-force) in a vac- 
uum region of the spacetime, the resulting equations of 
motion are ma 11 = f^ ct , where m is the particle's mass, a 11 
its covariant acceleration, and [LRR Sec. 2.5.6, Eq. (481)] 



This depends on the particle's entire history from r' — 
— oo to t' = +oo. Because the retarded and advanced 
fields are equally singular on the world line, the eval- 
uation of the radiation-reaction force does not require 
regularization. As was mentioned in Sec. II of the main 
text, the total work done by the radiation-reaction force 
equals the total energy radiated by the charged particle 



The expressions (|A3|) . ljA5|) . and 1A6(I for the retarded, 
advanced, and radiation-reaction forces are valid for any 
charge that moves freely in a vacuum region of an arbi- 
trary spacetime. Let us now specialize to the case of a 
particle moving slowly in the weakly curved spacetime of 
a (noncompact) star of mass M. The steps required to 
compute the retarded Green's function for such a space- 
time, and to evaluate the integral of Eq. (|A3|I . are de- 
tailed in PP [2fj. They lead to [PP Eq. (5.23)] 



, q 2 M „ 
/ret = 5-r 



2 2 dg 
3 Q dV 



(A7) 



which is also Eq. (j3Jl of Sec. IV, where the notation is in- 
troduced. Going through the same steps, starting instead 
with Eq. (|A5|) . gives 



/£ t (r) = 2q 2 u v f ^Gj x ,(z(r),z(r'))u y (r')dr' 

(A3) 

is the true, retarded, electromagnetic self-force. Notice 
that the integration extends over the particle's past his- 
tory, and that it is cut short at t' = t~ = t — + to avoid 
the singular behavior of the retarded Green's function at 
coincidence. This limiting procedure was first derived by 
DeWitt and Brehme [2{|, and it is a natural byproduct 
of the Detweiler- Whiting regularization method 0. 

The advanced version of the self-force is obtained by 
starting with the advanced solution for the potential, 



^adv( x ) 



G 



adv 



(x, z{t))u^{t) dr, 



(A4) 



and going through the same calculational steps as de- 
scribed above. The advanced electromagnetic field 
F"f (x) also is singular on the world line, but it is reg- 
ularized by the same singular field F^(x) as the re- 
tarded field. The end result for the equations of motion 
is ma^ = /^ dv , with 

/OO 
+ ^Gjl y (z{T),z{r'))u^\r')dT' 

(A5) 

being the advanced version of the electromagnetic self- 
force. Notice that the integration now extends over 
the particle's future history, and that it is cut short at 
t' = t + = t + + to avoid the singular behavior of the 
advanced Green's function at coincidence. 

Mine's radiation-reaction force is defined by [l^ 



ret 



ft 



adv 



(A6) 



adv 



2 2 dg 

3 9 df 



(A8) 



Not surprisingly, the advanced self-force is obtained from 
the retarded self-force by time reversal, t — > — t. Finally, 
from Eqs. (|A6|I - HA8|I we obtain Mino's radiation-reaction 
force, 



f - 2 
Jrr " 3 9 db' 



(A9) 



This is the familiar result that follows from the 
Abrahams-Lorentz-Dirac equation in flat spacetime (see, 
for example, Ch. 16 of Ref. [3^1- 



APPENDIX B: ORBITAL EVOLUTION UNDER 
THE SELF-FORCE 

In this Appendix we provide a derivation of Eqs. JSJ) 
. We rely on standard analytical methods of celestial 
mechanics as described, for example, in Ref. [3|J. 

We work in the x-y plane and introduce the polar co- 
ordinates (r, 4>) defined by x = r cos (f> and y — r sin <fi. 
We use the basis (r, (/>), in which r is a unit vector point- 
ing in the direction of increasing r, and </> is a unit vec- 
tor pointing in the direction of increasing cj>. The par- 
ticle's velocity vector, for example, can be decomposed 
as v = ff + {r<j})<p, where a dot indicates differentiation 
with respect to t. 

The self-force of Eq. Q admits the decomposition 



/self = m(Rr + Sep), 



(Bl) 
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with 



and 



R — ( A C ^T + ArrTr - ^ 

m \ r 3 



(see, for example, Ch. 12 of Ref. |23). In this method 
we average the right-hand side of each equation over a 
(B2) complete orbital period, treating (p,e, uS) as constants 
within the integral. If, for example, I(p,e,ui;v) denotes 
the time derivative of orbital element /, then we calculate 



q 2 M 



-A ri 



3r 2 



(B3) 



We treat q 2 /m as a small parameter and the self- force 
as a small perturbing force that produces only a slight 
deviation with respect to Keplerian motion. 

The charge's Keplerian orbit is described by Eq. 10}, 
which we rewrite as 



P 



1 



e cos v 



(B4) 



in terms of the true anomaly v; the longitude is then 
given by (f> = v + uj. For Keplerian motion the orbital 
elements (p, e,u>) are constant, and t(v) is determined by 
integrating the differential equation 



' M I, : - 

— 1 1 + e cos v ) 

pi \ 



Combining Eqs. (|B4|) and (|B5|) produces 



(B5) 



(B6) 



The orbital period is the time required by the particle to 
go from periapsis (v = 0) to apoapsis (v = n) and then 
back to periapsis (v = 27r) . It is given by 




P = 2tt* 



P 



(l-e 2 ) 3 M' 



(B7) 



During this interval 4> increases from u> to lo + 2tt. 

The electromagnetic self-force perturbs the particle's 
motion away from its Keplerian orbit. In the method of 
osculating orbital elements the perturbed motion is still 
described exactly by Eqs. (|B4j) - l)B6|l . but the elements 
(p,e,uj) acquire a time dependence. The driving equa- 
tions are (see, for example, Ch. 11 of Ref. [3^1 




S 



M 1 + e cos v ' 
R sin v + S 



-Rcosv + S 



2 cosv + e(l + cos 2 v) 
1 + e cos v 
(2 + ecosv) sinv 



1 



e cos v 



(B8) 
,(B9) 
,(B10) 



where R and S are listed in Eqs. (|B2|) and (|B3() . respec- 
tively. 

For the purpose of obtaining the secular behavior of the 
orbital elements, it is sufficient to integrate Eqs. (|B8|I - 
(|B10f) approximately, employing the method of averaging 



I(p,e,uj;v)dt = 



I(P, e,u;v) 



dv, 
(Bll) 

using Eq. I|B5|) for v and keeping (p, e, ui) fixed while in- 
tegrating. The result discards the rapid oscillations in 
I(t) which do not accumulate over a large number of or- 
bits. It keeps, however, the secular drift that eventually 
produces a large cumulative change. 

Substitution of Eqs. HB8(1 - 1|B10(I into the integral of 
Eq. I|B11(I leads to the results displayed in Eqs. JSJ 0. 
We recall that only the radiation-reaction component of 
the self-force is responsible for the secular evolution of p 
and e, while only the conservative component is respon- 
sible for the secular change in u>. 

We note that in the case of a purely conservative self- 
force [obtained by setting A rr = in Eq. (J3J] the relation 
between <f> and v can be obtained exactly: <p(v) = [1 — 
q 2 / {mp)] 1 ^ 2 v. This exact expression leads to a modified 
form of Eq. 0, 



2\3/2 




/(*), 



where x = q 2 j(mp) and 



/(*) = 



x/2 



= 1 



x 
1 



(B12) 



(B13) 



When q 2 <^ mp we have that f(x) ~ 1 and Eq. Q is a 
very good approximation to Eq. I|B12|I . 

In Appendix C we will need explicit expressions for 
(p)(v), (e)(v), and (ui)(v). To obtain these we first aver- 
age I' = dl/dv = I/v over a complete orbital period, to 
get (/') ee (2TT)' 1 I' dv = (P/2tt)(7). This yields 



Aq 2 IM 
-A rr — — -i / — , 
3 m 



(P) 



(e) = -A rr — e 

m \i p- 




2 rap 



Next we integrate, and obtain 

(P) = Po(l - KrV/v ) 2/3 

(e) = e (l - \ n -v/v ) 1/2 

3 fpa_ 
4V M 



{a;} 



1 - (1 - X c v/v Q ) 



1/3 



(B14) 

(B15) 
(B16) 

(B17) 
(B18) 

(B19) 
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where 




(B20) 



We assume that the orbital evolution starts at v = 
with values (p — po,e = e$,uj = 0) for the orbital el- 
ements. The evolution ends at v = v when both (p) 
and (e) are zero; the final value of the pcriapsis shift is 
(lu){v ) = -(3/4)y/p /M, and this expression confirms 
the estimate of Eq. itTUI) . 



APPENDIX C: ELECTROMAGNETIC WAVE 

We calculate the electromagnetic wave (or more pre- 
cisely, the wave-zone magnetic field) produced by a 
charged particle subjected to the electromagnetic self- 
force of Eq. As was explained in Appendix B, the 
motion is at all times described by r(v) given by Eq. (|B4|1 . 
4>(v) = v + lo, and t(v) is obtained by numerically inte- 
grating Eq. I|B5I) . The orbital elements (p, e,uj) evolve 
according to Eqs. l|BT7jl - l|B20)l . 

In the electric-dipole approximation the wave-zone 
magnetic field is given by (see, for example, Sec. 11.1.4 
of Ref. H|) 



B 



47T cd 



(CI) 



where h = [sin a cos /3, sin a sin /3, cos a] is a unit vector 
which points from the source to the detector, and p is 
the source's electric dipole moment expressed in terms 
of retarded time t — d/c; the detector is at a distance 
d from the source and c is the speed of light. For an 
orbiting charge we have p — qa, where a is the parti- 
cle's acceleration. For our purposes it is sufficient to use 
the approximation a = g = — (M/r 2 )r and Eq. (|C1|) 
becomes 



B 



po qM n x r 
47r cd r 2 



(C2) 



With f = [cos 4>, sin 0, 0] we have 



n x r = [— cos a sin <j>, cos a cos cj>, sin a s'm(4> — 0)], (C3) 

and each component of the magnetic field can be com- 
puted straightforwardly. 

To produce the plots of Fig. 1 we chose the viewing 
angles (a — 7r/4, (3 — 0) and selected the y component of 
the magnetic field. To integrate the equations of motion 
we set po = 20M, e = 0.5, and q 2 /(mp ) = 0.002. The 
magnetic field was rescaled by an arbitrary numerical 
factor to obtain values of order unity. To produce the 
true wave (represented as a blue/dark curve) we set A rr = 
A c = 1 in Eqs. lB17(l - ilB19(l . To produce the adiabatic 
approximation to the wave (represented as a red/light 
curve) we selected A rr = 1 and A c = 0. 
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